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ABSTRACT 


Based  upon  the  concept  of  Galerkin's  approximate  method  for  solving 
eigenvalue  problems,  a  new  scheme  of  numerical  treatment  is  proposed  for 
a  class  of  nonconservative  (circulatory)  elastic  stability  problems.  This 
is  accomplished  by  considering,  together  with  the  original  system,  a 
second  system  which  is  obtained  by  introducing  an  adjoint  to  the  circulatory 
force  field.  The  resulting  problem  is  shown  to  be  self-adjoint,  with  eigen¬ 
functions  which  possess  the  property  of  reducing  the  original  problem  to  a 
simple  integral  equation.  This  Integral  equation  may  be  solved  by  quad¬ 
rature  and  an  estimation  of  error  is  also  possible.  The  proposed  method 
is  especially  suitable  for  direct  evaluation  on  a  digital  computer  and 
does  not  involve  tedious  integrations  of  functions  encountered  in  the 
commonly  adopted  application  of  the  Galerkln  method. 


INTRODUCTION 


Galerkin's  method  has  proved  to  be  a  powerful  technique  for  constructing 
approximate  solutions  of  boundary  value  problems  of  engineering  science. 

Following  Kikhlin  [1],  a  solution  is  assumed  in  the  form  of  a  known  set  of 
functions  with  unknown  generalized  co-ordinates  or  constants.  The  series, 
when  substituted  in  the  differential  equation,  gives  rise  to  an  error  or  residue 
which  is  required  to  be  orthogonal  to  the  members  of  the  set.  This  results  in 
a  set  of  equations  for  determining  the  unknowns.  The  dev.  lopment,  In  recent 
years,  of  several  closely  related  methods,  such  as  the  least  square  technique, 
the  collocation  method,  and  the  method  of  moments,  indicates  the  success  and  a 
general  sense  of  confidence  in  the  application  of  the  concept  of  Galerkin. 

This  concept  was  originally  introduced  for  solving  self-adjoint  boundary 
value  problems  for  which  a  general  proof  of  convergence  hus  been  supplied  by 
Mikhlin  [1].  Bolotin  [2]  has,  however,  suggested  that  by  taking  the  eigen¬ 
functions  of  f i “ft  vibration  as  the  coordinate  functions,  nonself-adjoint  problems 
of  elastic  stability  with  follower-type  loads  may  also  be  approximately  solved 
by  the  Galerkin  method.  In  this  way  he  has  successfully  treated  a  number  of  non¬ 
conservative  stability  problems.  But  it  appears  that  more  effective  generaliza¬ 
tions  of  the  Galerkin  method  in  this  area  of  nonconservative  elastic  stability 
are  possible,  and,  in  particular,  studies  which  dwell  on  the  question  of  convergence, 
estimation  of  error,  and  direct  efficient  use  of  digital  computers  are  of  interest. 

The  present  study  discusses  an  extension  of  the  Galerkin  method  for  direct 
numerical  treatment  of  a  class  of  nonconservative  stability  problems.  For  this 
purpose  we  consider  an  additional  system  by  combining  the  force  field  of  the 
adjoint  to  the  original  problem.  The  authors  have  recently  proposed  variational 
methods  in  which  the  concept  of  adjoint  systems  was  employed  [3,^].  We  prove 
first  that  for  a  certain  class  of  systems  the  resulting  boundary  value  problem 
is  self-adjoint.  Therefore,  although  both  fields  of  forces  are  separately 
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nonconservative,  the  combination  of  these  is  a  conservative  system  of  forces. 
This  feature  was  first  noted  by  Nemat-Nasser  and  Herrmann  [5'i  who  combined  the 
force  fields  of  Beck's  and  Reut's  models  f2j.  This  seif-adjoint  system  bears 
close  relationships  to  the  original  system,  such  that  in  the  Galerkin  method 
the  coordinate  functions  may  be  selected  as  being  the  eigenfunctions  of  this 
resulting  system.  By  means  of  this  set  of  eigenfunctions,  it  is  possible  to 
transform  the  boundary  value  problem  into  an  integral  equation  which  may  be 
solved  numerically  by  quadrature. 

The  essence  of  the  proposed  extension  of  the  Galerkin  method  for  the 
class  of  nonself-adjoint  system  selected  for  study  may  be  described  as  fc.(low= 

In  an  extended  sense,  the  Galerkin  method  requires  the  solution  of 

<Lu  ,  v>  ■  <1  ,  v> 

for  the  problem  Lu  »  f  with  appropriate  boundary  conditions,  where  the  un¬ 
known  functions  u  is  to  be  approximated  by  the  coordinate  functions 
If  we  select  v  such  that  the  above  may  be  reduced  to 

■A’ 

<L  v  ,  u>  «  <f  ,  v> 

then  this  integral  equation  mo,;  i?  solved  by  means  of  3  quadrature  whose  con¬ 
vergence  property  is  known  and  certain  error  estimations  may  also  be  given.  The 

A: 

nature  of  L  ,  the  properties  of  v  and  certain  quadrature  schemes  are 
discussed  below. 


Please  note  the  difference  between:  L  (capital  t) 

and 

L  (script  t) 


AJ1ASS,  COONCQNSERVATIVE  STABIL1JY.  PROBLEMS 


We  consider  the  following  form  of  an  ordinary,  linear,  differential 
equat ion: 


Lu  +  BL  u  *  Xu 


0) 


where 


N  ,n 

L  =»  7  au  (x)  - — 

i  N-n  ,  n 

n«l  dx 


(2) 


l 


K  <  N 


(3) 


n=*l 


dx 


In  the  above  u  denotes  a  function  of  a  real  variable  x  for  a  <_  x  b  , 
c»n  are  continuous  functions  of  x  whose  M-n  derivatives  with  respect  to  x 
exist  ar.d  are  continuous,  and  6n  are  certain  constants.  B  is  a  parameter 
representing  the  magnitude  of  the  follower  load.  Further,  aQ  does  not  vanish 
at  any  point  of  the  closed  interval  (a  ,  b) . 

Associated  with  (1)  we  consider  N  linear,  homogeneous,  boundary  con¬ 
ditions  in  u(a)  ,  u' (a)  ,  .  ,  u  ^  (a)  ,  u(b)  ,  .  ,  u^N  ^  (b)  ,  as 

given  by 


L.u 

J 


,n 

d  u 


r  a  u 

njn  dxn 


n«*0 


j  -  1  ,  2  ,  ...  ,  N 


(4) 


njn  are  quantities  characterizing  certain  properties,  such  as  stiffness  or  in¬ 
ertia  at  the  end  points  (a  ,  b) .  For  future  use  let  us  define  N  additional 


forms  L^+jU  ,  . L^u  in  u‘(a)  and  u 1  ( b)  ,  i  =0  ,  1  . .  N-l  , 


so  that  LjU  ,  L2u  , 


,  l2(s|U  are  linearly  independent. 


1» 

In  several  nonconservative  structural  stability  problems  it  turns  out 
that  the  follower  load  influences  the  governing  equations  of  motion  by  the 
presence  of  the  operator  L in  (l)  so  that  in  the  absence  of  such  loadings 
the  free  vibration  of  the  system  is  expressed  by  the  same  operator 

Lv  ■  uv  (5) 


with  the  same  boundary  conditions 

IjU  *  0  ,  j  »  1  ,  2 .  ,  N  (6) 

The  eigenvalue  problem  governed  by  ( 5“ 6)  Is  known  to  be  self-adjoint  and  thus, 
the  presence  of  the  operator  Lj  in  (1)  destroys  this  property.  It  is  well  to 
emphasize  that  the  class  of  systems  selected  for  study  is  such  that  the  follower 
load  does  not  affect  the  boundary  conditions.  It  is  further  noted  that  this 
class  admits  the  operator  L1  with  constant  coefficients.  In  a  later  section 
we  will  Investigate  the  problem  by  rela  Ing  this  restriction  and  it  will  be 
shown  that  several  systems  with  Lx  having  variable  coefficients  also  possess 
the  same  property. 

From  (1)  we  obtain,  after  integrating  by  parts  t6], 

<u*  ,  (L  +  8L1)u>  -  <u  ,  (1*  +  j )u*> 

■  [P(u  ,  U*)]b  +  8 [P  (u  ,  u*)lb 

a  a 


b 

<u  ,  v>  =  /  uvdx 
a 


where 


5 


denotes  spatial  average  and 


.*  *  ,  ,.N 

i  u  «  (-1; 


.N,  \ 

d  (anu  ) 


.N-l,  *\ 

0“  '  .  ,  .,N-1  d  (ctiU  ' 


dx 


N 


+  H)n_1 


dx 


N-l 


+  . . 


d(aN_1u") 

dx 


(8) 


lV«(-1)KS  (-D^b  dK”^u* 

L>u  (,)  '•*,«  +  ( *'  8>^  + 


-8 


K-l 


du 

dx 


,  K  <  N  (9) 


and  [P (u  ,  u  )]b  and  [P,(u  ,  u  )]b  are  bilinear  forms  in  u(a)  ,  u 1  (a)  , 
.  ,  u(N-,)(a)  ,  u(b)  ,  .  ,  u{N_,)(b)  ,  u*(a) . u*(N_1)'b)  . 


By  insisting  on  the  property  of  system  (5-6)  to  be  self-adjoint,  it  may 
be  stated  thfit  P  must  be  expressed  in  the  following  form: 


[P(u  ,  u*)]b  -  l  LuL 

a  i  n 


2N 


n™l 


2N+L-n 


u 


(10) 


and  further  that 


(ID 


An  examination  of  (8)  reveals  that  a  necessary  condition  for  equation  (U)  to  be 
valid  is  that  N  must  be  an  even  integer,  in  other  wods,  the  order  of  the 
differential  equation  must  be  even. 

It  is  now  possible  to  properly  select  2K  linearly  independent  forms  M,  , 


M2  ,  ....  ,  in  u*(a)  ,  u*5  (a) ,  -  ,  u  (a)  ,  u’(b)  , 

so  that  P  may  be  written  in  the  following  form 


u*(K“,J(b) 


2K 


[P,  (u  ,  u'')]b  -  l  L  u  M_  .  u"  ,  K  < 

1  ’  a  *\  n  2K+l-n  ’  — 

n-l 


(12) 


Therefore,  an  adjoint  to  system  (l-4*)  may  be  defined  by  the  following  eigen¬ 
value  problem 


6 


*****  .... 

Lu  +  BL u  «Xu  (13) 

with  the  following  boundary  conditions 

LjU*  +  BMjU  ■  0  ,  j  ■  1  ,  2  ,  ,  N  (14) 

Note  that  M^+j  ■  **  0  . 

It  has  been  shown  in  past  studies  [3,4],  that  the  two  sets  of  eigenvalues 
{ X | }  and  {X.*}  are  identical.  Note  that  the  system  (1-4)  will  be  self- 

/V  & 

adjoint  only  if  MjU  =  0  and  “  Lj  ,  which  is  not  the  case  for  noncon- 

servative  stability  problems. 

LSZLFsmom  SYSTEM 

Keeping  in  mind  the  developments  oF  the  previous  section,  one  may  hope  to 
recover  another  system  which  is  self-adjoint  and  whose  eigenvalues  have  close 
relationships  with  those  of  the  original  problem.  L,  represents  the  field  of 

.L 

the  follower  load  whose  adjoint  force  field  is  represented  by  L j  .  Generally, 
this  concept  of  the  adjoint  force  field  in  nondissipative,  nonconservative 
systems  is  quite  abstract  and  its  physical  significance  mny  be  traced  only  in 
simple  problems.  Indeed,  such  a  relationship  between  Beck's  and  Reut's  models 
was  discovered  with  some  surprise  by  Nemat-Nasser  and  Herrmenn  [5]  who  also 
pointed  out  that  a  combination  of  the  two  nonconservative  forces  results  in  a 
conservative  system.  In  this  section  we  examine  a  similar  property  mathemati¬ 
cally  for  the  group  of  problems  under  consideration. 

We  define  an  eigenvalue  problem  by  the  following  differential  equation 

a. 

Lv  +  BljV  +  BL.  v  =  u)V 


05) 
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with  the  boundary  conditions 


Uv  +  gMjV  =  0  ,  j  »  1  ,  2  , . N  (16) 

and  prove  that  it  is  a  self-adjoint  system. 

After  integrating  by  parts,  one  obtains 

<v  ,  (L  +  BL  +  BL  )v>  *<v  ,  (L  +  bL  +  BL  )v* >  +  [P (v  ,  v  )]^ 

11  1  1  9 

+  B[P1(v,vXb*6[P2(v.v*)]ab  <17) 


The  bilinear  concomitants  Pj  and  P2  may  be  written  explicitly  as 
Pj(v  ,  v*)  -  vtS^v*  -  SK.2^+  .  +  (-')K''Bo 

+  6k-3  —  + . +  (-l)K-26o 

dx  dx 

dK''v  „  * 

+ . +  -  K—  V 

dx 

and 

P2(v  .  V4)  -  v[-6k.,v*  -  Bk.2  . * 

*  7  [«K.2»‘  +  V-3  f- * . M-l)2K-\  gf] 

„  .K-l 

+ . *  <-')K7tptM‘ 


(18) 


(13) 
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Before  we  proceed  further,  it  may  be  noted  that  if  the  operator  consists 
of  only  odd  derivatives,  then 


L1  +  L1 


(20) 


and 


P^v  ,  v*)  +  P2(v  ,  v*)  -  0  (21) 

In  this  special  case  a  proper  eigenvalue  problem  would  be  defined  by  Lv  ■  wv 

with  the  boundry  conditions  LjV  ■  0  ,  j  ■  1  ,  2  . . N  ,  which,  incidentally, 

coincides  with  the  equation  of  motion  of  free  vibration  and  is  obviously  self- 
adjoint.  It  may  be  further  observed  that  in  a  general  situation  (Lj  +  L1  ) 
together  will  cancel  all  odd  derivatives  and  will  consist  of  terms  with  only 
even  derivatives.  Similarly,  {?x  +  P2)  will  also  cancel  terms  with  coefficients 
»  •••  and,  therefore,  in  the  following  prooP  we  will  insist  that 
L  consists  of.  only  even  derivatives  and  K  is  an  even  integer.  Th"e,  a  close 
examination  of  ( 1 8- i^)  reveals  that 


P,  -  P„  (22) 

1  2 

and  (17)  reduces  to 

.  >v  .  * 

<v  ,  (/  +  pLj  +  BL1  )v>  «  <v  ,  (L  +  tiL  +  BL^v  > 

+  [P(v  ,  v*)]*  +  28 [P  (v  ,  v*)]*  (23) 

d  a 

The  bracketed  terms  on  the  r;^ht  side  of  (23)  may  now  be  expressed  as 
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*  k  ,  ,,  211 

IMv  ,  v  )]a  ♦  28[Pj(v  ,  v")]b  -  I  (L„  *  8Hn)v(L2Ht,.a  ♦  BM2N+,.„)»'  (24) 

n*  i 

and  (23)  yields 

<v*  ,  (L  +  bL  +  6L1*)v>  -  <v  ,  {L  +  BL1  +  B L*)v'> 

2N 

■  l  <L„  +  BM„>''<L2N+l-n  +  SWn)v'  (2f) 

na  i 

Thus  if  (16)  is  obeyed  by  v  and  v  ,  the  right  hand  side  of  (25)  is  zero. 
Therefore,  the  theorem  is  proved. 

When  operator  L consists  of  both  odd  and  even  derivatives,  a  proper 
selection  of  a  self-adjoint  system  would  consist  of  the  same  equations  (15)  and 

(16)  with  the  restrictions  that  terms  involving  Bk_j  ,  ,  .  are 

identically  zero.  This  follows  in  a  natural  way  from  the  above  proof. 


in  the  previous  section  it  was  possible  to  study  the  properties  of  the 

bilinear  concomitants  P,  and  P  in  explicit  forms  which  facilitated  the 

1  2 

proof.  For  a  more  ge  jral  operator  Lj  having  variable  coefficients,  it 
will  be  seen  now  that  an  indirect  approach  must  be  adopted.  A  wide  group  of 
problems  in  structural  stability  with  distributed  follower  loadings  gives  rise 
to  the  study  of  eigenvalue  problems  in  which  the  operator  ^  involves  var¬ 
iable  coefficients.  It  has  been  experienced  that  distributed  follower  loading 
is  not  the  only  source  of  the  origin  of  variable  coefficients,  but  several 
other  problems  such  as  the  instability  in  bending-torsion  of  a  rectangular  bar 
acted  upon  by  a  transverse  follower  force,  also  have  a  similar  nature.  Thus, 
it  seems  desirable  to  extend  the  results  cf  the  previous  section  to  encompass 
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a  wider  group.  Accordingly,  we  consider  the  following  form  of  l  along  with 
(1)  and  (4) 


K  n 

Ll  "  i.  eK-n(x>  TTT  •  KiN  (2« 

n»1  dx 

The  8n(x)  's  are  continuous  functions  of  x  whose  K-n  derivatives  with 
respect  to  x  exist  and  are  continuous. 

We  note  the  following  relation: 


<u 


Llu> 


JL  JL 
/T 


<U  ,  Lj  u  >  +  [P  (u  ,  u  )]. 


(27) 


where 


*  * 


Lj  u  -  (-1) 


,K-1 to 

K  “  '"0“  '  .  ,  , XK-1  d  lBlU  } 


d*(Bnul 


dx 


+  (-1) 


d(BK-lu") 


dx 


(23) 


dx 


and  P3  Is  a  bilinear  form  in  u(a)  ,  u‘ (a)  .  u^K_,^(a)  ,  u(b)  ,  . 

^  (b)  ,  u  (a)  . . u  (a)  . .  u  ^  ^  (b)  .  Analogous  to  the  previous 

case,  P3  may  be  expressed  as 


(P3(u  *  UX  3  l  LnuM2K+1-nl 
n=“l 


2K 


(29) 


where  Mj  ,  M2  .  MjK  are  linearly  independent  forms  in  u  (a)  and 


*i  . 

u  Cb)  ,i—0,l,  . K  .  On  the  other  hand,  we  have 


*  .  ?V  ,  *  r  /  * ,  •,  b 

<u  ,  L1  u>  «  <u  ,  LjU  >  +  [P4(u  ,  u  )la 


(30) 


The  bilinear  concomitants  P,  and  Pv  have  the  following  forms  in  u1  and  u  1 
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,  *x  r  *  d(6K-2U  } 

P  fu  ,  u  )  -  u[B„  .u - — - + 

3  *■'  dx 


,  dK*'  <8 /) 

♦  <-'>  -ttfH 

dx 

jK-2/„  *, 


du  *  d(B,,  )  K_ 2  d  (S0u  ) 

+  —  [e..,u  -  ♦ . ♦  (-:)K  2  — A— 


jK-1 

d  u  * 
+  “kT  60U 

dx 


*  *  J(BK-2U)  KdK,(60U^ 

P  (u  ,  u  )  «  u  [-&,  .u  + - + . +  (-1)K - jrn — ] 

4  K  1  dx  dxK  1 

du*  d^6k  k  i  dK  2^enu  ^ 

♦  —  t-6K-2u  +  ~~~~  + . *  (-'>  —id—' 1 

dx  dx  dx 

*  . 


Thus  we  see  that  (32)  may  be  obtained  from  (31)  after  replacing  u  by  u  and 
u  by  -u  ,  and  vice  versa.  Consequently,  a  counter  part  of  (29)  may  be 
written  for  P,  bv 


r  ,  *.,b  r  * 

[p4<u  •  u  )!a  ■  \  VL2K+I-4U 
n-1 


We  now  show  that  the  eigenvalue  problem  given  by  the  differential  equation 


Lv  +  BLjV  +  SL 1  v  ■  ojv 


with  the  boundary  conditions 


L .v  +  BM.v  =  0  , 

J  J 


j  -  1  ,  2  , 


is  self-adjoint. 
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As  we  have  discussed,  one  obtains 

<v"  ,  (L  +  BLj  +  8L1*) v>  =  <v  ,  (L  +  +  [P(v  ,  v*)]^ 

+  8 [P 3 (v  ,  «*)]J  +  BfPjv  ,  v*)]J  (36) 

which  yields,  after  using  (29)  and  (33) 

<v‘  ,  (L  •{-  BL1  +  Bl-1  )v>  -  <v  ,  (L  +  Blj  +  SLj  )v*> 

2N 

'  j,  (L"  +  6Hn)v(L2N*l-n  +  6fWn>v'  <”> 

and  if  (35)  is  obeyed,  the  right  hand  side  of  (37)  vanishes.  This  concludes 
the  proof. 

RELATIONSHIPS 

In  view  of  the  fact  that  the  system  governed  by  (3^-35)  is  self-adjoint, 
it  is  capable  of  yielding  accurate  solutions  based  upon  the  usual  Galerkin 
procedure  or  other  approximation  methods.  Therefore,  it  is  of  interest  to 
investigate  if  this  system  has  any  '•elationship  with  the  original  problem. 

Indeed,  It  will  be  seen  in  the  next  section  that,  based  upon  certain  relation¬ 
ships  discussed  in  the  following,  an  extension  of  the  Galerkin  method  for  the 
solution  of  the  original  nonconservative  system  may  be  advanced  which  has 
various  desirable  features. 

Let  us  consider  the  following  boundary  v=ilue  problems:; 

Lu  +  =  f  (38a) 


L.u  =  0 
J 


j  =  1  ,  2  .  ,  N 


(38b) 
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and 

Lv  +  Bl^v  +  U^v  •  f  (39a) 

LjV  +  $H.v  -  0  ,  j  »  1  ,  2 . N  (39b) 

where  f  Is  an  arbitrary,  continuous  function  of  the  variable  x  .  Multiply 
(38a)  by  v  and  integrate  in  the  domain  of  x  to  obtain 

<v  ,  (L  +  3l-1)u>  -  <u  ,  (L  +  BL>>  *  <v  ,  f>  (^0) 

which  yields  after  using  (39) 


<U  ,  *  <u  ,  f>  -  <v  ,  f>  (4l) 

Similarly,  starting  with  (39a)  and  using  (38)  we  obtain 

<v  ,  BLj  u>  ■  <u  ,  f>  -  <v  ,  f>  (42) 

Thus  by  comparing  (41)  and  (42)  we  have 

<u  ,  LjV>  -  <v  ,  41*u>  (43) 

which  is  reminiscent  of  the  Betti-Rayl eigh  reciprocity  relation  in  classical 
solid  mechanics.  Note  that  the  operator  L  gives  the  field  of  nonconservative 
loading  and  a  spatial  average  of  the  type  <u  ,  L^v>  represents  work  done  by 
LjV  through  displacement  u  . 
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It  may  be  pointed  out  here  that  a  more  direct  relationship  exists  with 
the  adjoint  boundary  value  problem.  Along  with  (38)  consider  the  following 

Lu"  +  SL^u*  -  f  (44a) 

L.u*  +  8H.U*  -  0  ,  j  -  1  ,  2 .  ,  N  (44b) 

Multiplying  (38a)  by  u  and  integrating  in  the  domain  of  x  ,  we  obtain 

<u  ,  (L  +  SL1)u>  ■  <u  ,  (L  +  pt1*)u*>  *  <u*  ,  f>  (45) 

If  we  now  use  (44a),  (45)  yields 

<u  ,  f>  *  <u*  ,  f>  (46) 

which  states:  The  spatial  averages  of  the  products  of  the  external  loading  with 
the  displacements  of  a  nonconservative  system  and  its  adjoints  are  the  same. 


In  the  theory  of  stability  wide  use  is  made  of  approximate  methods  to  de¬ 
termine  the  regions  of  stability  of  a  particular  problem.  Accurate  determina¬ 


tion  of  critical  parameters  which  separate  stable  from  unstable  regions  of  a 


conservative  system  may  be  facilitated  by  the  use  of  a  particular  technique 
available,  (e.g.  Methods  of  Ritz,  Rayleigh,  Galerkin,  etc.)  all  of  which  have 
their  origins  in  variational  principles  [1]. 


Several  extensions  of  these  methods  have  been  suggested  for  nonconserva 


tive  systems  which,  however,  are  only  heuristicsl ly  justifiable  and  lack  proofs 
of  convergence  and  estimations  of  error  involved.  Most  authors  analyze 
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nonconservative  stability  problems  in  two  parts.  First,  a  related  simpler 
eigenvalue  problem,  (usually  governing  equations  of  free  vibration  which  con¬ 
stitutes  a  self-adjoint  '/stem,  see  Bolotin  [2]), which  yields  a  complete  set 
of  trial  functions,  is  solved.  Then,  with  the  help  of  this  set  of  trial 
functions,  an  application  of  the  approximate  procedure  adopted  reduces  the 
problem  to  determining  the  roots  of  a  certain  characteristic  determinant.  How¬ 
ever,  since  no  firm  guide  line  in  the  selection  of  the  trial  functions  existed, 
the  establishment  of  convergence  and  an  estimation  of  the  error  involved  is  not 
readily  possible. 

In  the  following  we  suggest  an  extension  of  the  Galerkin  method  of  solving 
the  general  eigenvalue  problem 


Lu  +  BLjU  ■  Xu  (46a ) 

LjU  ■  0  ,  j  -  1  ,  2 .  ,  N  (46b) 


using  the  eigenfunctions  of  the  self-adjoint  system  defined  by 


A 

Lv.  +  BL.v.  +  $L  v.  -  w.v.  (47a) 

J  1  J  i  J  J  J 

L  v,  +  BM  v.  ■  0  ,  m  ■  1  ,  2 .  ,  N  (47b) 

m  j  m  j 

We  denote  by  {to.}  and  { v j }  the  sets  of  eigenvalues  and  eigenfunctions,  re¬ 
spectively,  of  (47)  and  consider  the  following  integral 


<v . 
J 


(Lu  +  BLjU)>  *  X < v j  ,  u> 


which,  after  integration  by  parts  and  in  view  of  (1*7),  reduces  to 
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<u  ,  BL^v.*  -  (ok  -  X)<u  ,  v.>  (48) 

As  is  evident  from  (48),  we  Vive  converted  the  differential  equation  to 
the  simplest  kind  of  an  integral  relation.  In  fact,  this  integral  equation 
necessitates  determining  the  area  under  a  curve  given  by  he  unknown  function 
u  which  is  multiplied  by  certain  known  functions  Involving  Vj  and  their 
derivatives.  Thus,  in  accordance  with  practice  in  numerical  analysis,  the  inte¬ 
gration  may  be  replaced  by  the  summation  of  a  finite  series  in  terms  of  the 
ordinates  of  the  function  u  at  certain  intervals.  There  are  several  well 
known  rules  of  carrying  out  this  type  of  numerical  integration  for  which  an 
estimate  of  the  round-off  error  may  be  obtained.  In  order  to  illustrate  the 
procedure  of  the  proposed  method,  we  obtain  in  the  following  a  general  quadrature 
relation  based  on  Newton's  interpolation  formula.  Let  u  take  the  values  uQ  , 
u1 .  ,  uN  for  the  equidistant  points  xQ  ,  . . xN  of  the  inde¬ 

pendent  variable  x  with  xQ  ■  a  and  x^  ■  b  ,  the  interval  being  h  .  In 
Newton's  formula  we  have 

x  ■  x  +  hn 

0 

from  which  we  get 

dx  -  hdn 


+ 
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where 


Au  -  u  -  u. 
0  1  o 


A2u0  -  AUj  -  AUQ  -  u2  -  2Uj  +  uQ 


3i'o  ■  a2u!  "  A2uq  -  u3  -  3u2  +  3u3  -  uc 


etc. 


Substituting  (1*9)  in  (48)  and  carrying  out  the  integration  with  respect 
to  n  ,  the  following  expression  is  obtained 


N2  N3  N2  A2uq  N4  A3uc 

f*aojuG +  -  aij4uo  *  <r  ‘  r)a2j-2T +  <r  • 1)3  +  ""Ntt 


N5  3N4  1 1N3 


A4u, 


+  ( I — r +  —  ‘  3N2)a4j  tt  + . 1 

N2  N3  N2  A2u  N4  A3u. 

<“j  -  x)[Nb0ju0  +  T  b’jAu0  +  (r ‘F)b2j  TT+  V’  N3  +  N2)b3j  TT 


N5  3N4  1  IN3 


A4u, 


+  <r-~ +  — •  3N2)b.j-^+ . 1 


(50) 


where 


a  .  ■  Qi  v .  (x  +  hn) 
nj  i  J  O 


and 


b  .  »  v, (xn  +  hn) 
nj  j  o 

Keeping  in  mind  that  A  u  may  be  expressed  in  terms  of  u  .  u,  ,  .  ,  u, 

equation  (50)  may  be  arranged  in  the  following  form 
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,N+1  (N+l)m 

EN  *  '('n ~+  TT7 ’  ~ '  n<n  ‘  D...-(n  -  N)  ,  a  <  3  <  b  (53) 

Thus,  we  see  that  the  round-off  error  may  be  made  as  small  as  we  please  by 
making  the  Interval  smaller.  Estimates  of  the  error  Involved  In  the  methods  of 
Simpson,  Weddle,  Gauss,  etc.  a-e  discussed  In  [7]  and  [8], 

The  above  scheme  of  numerical  treatment  does  not  require  Integration  of  any 
kind,  by  contrast  to  the  commonly  employed  approximate  methods  and,  therefore,  Is 
suitable  for  direct  use  on  a  digital  computer.  It  may  be  mentioned  that  the  same 
scheme  may  be  found  more  convenient  also  for  self-adjoint  systems, in  which  case 
v  In  (**8)  may  be  any  complete  set  of  functions  satisfying  the  same  boundary 
conditions  as  those  for  u.  In  a  forthcoming  study,  numerical  results  of  several 
nonconservative  stability  problems  wtll  be  presented. 
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